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Nuclear Magnetic Resonance studies of Ga stabilized 5-Pu reveal detailed information about the 
local distortions surrounding the Ga impurities as well as provides information about the local spin 
fluctuations experienced by the Ga nuclei. The Ga NMR spectrum is inhomogeneously broadened 
by a distribution of local electric field gradients (EFGs), which indicates that the Ga experiences 
local distortions from cubic symmetry. The Knight shift and spin lattice relaxation rate indicate 
that the Ga is dominantly coupled to the Fermi surface via core polarization, and is inconsistent 
with magnetic order or low frequency spin correlations. 
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Introduction 

The investigavtion of the low temperature properties 
of plutonium and its compounds has experienced a re- 
naissance in recent years, and several important experi- 
ments have revealed unusual correlated electron behavior 
Q, 1^ 1^ . The 5f electrons in elemental plutonium are on 
the boundary between localized and itinerant behavior, 
so that slight perturbations in the Pu-Pu spacing can give 
rise to dramatic changes in the ground state character. 
In fact, Pu exhibits six allotropic forms in the solid state, 
with varying degrees of symmetry. The inability for band 
structure calculations to predict all six of these forms is 
testament to the complexity of the correlated electron 
behavior 0, |^. (5-Pu, which exists between 576K and 
720K, is the most technologically important form, and 
can be stabilized down to T = by doping with Al, Ga 
or In, which substitute for the Pu sites [g. In particular, 
Pui-xGax possesses the fee structure and physical prop- 
erties of (5-Pu for 0.020 ^ a; < 0.085. It is not obvious why 
the small amount of the secondary element will stabilize 
the electronic structure of the S phase, although clearly 
the slight perturbation reflects the itinerant-localization 
instability of the parent material. 

Several theories have emerged to try to explain the sta- 
bility of the S phase. One technological challenge of band 
calculations in Pu is the ability to predict local versus 
itinerant behavior for the five 5f electrons, because the 
theoretical approaches for these two extremes are funda- 
mentally different. These models typically predict either 
non-magnetic ground states, or ones where local Pu mo- 
ments are retained 0,0] • In principle, such predictions 
can be tested by various experimental techniques. 

Nuclear magnetic resonance (NMR) is ideal for prob- 
ing magnetic correlations and local structure. ^^^Pu has 
a nuclear spin I — i, so in principle it can serve as a di- 
rect probe of the magnetic correlations via the hyperfine 
interaction A between the nuclear (/) and the electronic 
(5) spins. However, for 4f and 5f nuclei, A is usually 
quite large (on the order of 1000 kOe/^s), so that fluc- 



tuations of of the electron spin S relax the nuclei so fast 
that their signal is rendered invisible. On the other hand, 
the hyperfine coupling to the nuclei of the secondary el- 
ement (Al, Ga or In) are typically one to three orders of 
magnitude smaller, so by measuring the secondary nuclei 
one can gain considerable insight into the spin dynamics 
of the system. 

One of the challenges facing band theorists calculating 
the electronic structure of S-Pu is the role of the sec- 
ondary atoms. Recent x-ray absorbtion fine structure 
(XAFS) studies suggest that Pu lattice distorts locally 
around the Ga sites 0. Clearly, these local distortions 
are caused by the electronic system, and should be cap- 
tured by band structure calculations. Therefore, a de- 
tailed experimental study of these local distortions puts 
important constraints on any theory. NMR is also ideal 
for probing these local distortions, and serves as a com- 
plement to the XAFS data. Nuclei with spin I > \ 
experience electric quadrupolar interactions between the 
quadrupolar moment Q of the nuclei and the local elec- 
tric field gradient (EFG), Vqq, where V is the crystal 
field potential. Ga [I — |) is therefore sensitive to the 
EFG at the impurity site. 

In this paper we present Knight shift {K), spin lat- 
tice relaxation rate [T^^), and linewidth (cr) data on 
the Ga nucleus in a polycrystal sample of Pui_a;Gaa; for 
X = 0.017 between 4K and lOOK. In the early 1970's, 
Fradin, Brodsky, and coworkers at Argonne performed 
the first NMR experiments on (5-Pu stabilized with Al 
0. They found no evidence for any magnetic transition 
or the formation of local moments. Our work on Ga sta- 
bilized (5-Pu provide the only new data since the pioneer- 
ing work of Fradin, and we find similar behavior. The 
Knight shift and T^^ data behave as in a typical con- 
ductor, and we find no evidence for either local moments 
at the Ga sites, or for strongly temperature dependent 
spin fiuctuations . On the other hand, we do not find 
that the EFG at the Ga site vanishes, as expected for a 
site with cubic symmetry. Rather, the Ga spectrum is 
inhomogeneously broadened by a distribution of EFGs. 
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FIG. 1: The ^^Ga NMR spectrum at 4K. The sohd vertical 
hne is the position of pure Ga metal, and the solid line through 
the data is a Gaussian fit as described in the text. 



FIG. 2: The Knight shift (•) and susceptibility (solid line) of 
the Ga versus temperature. The susceptibility data are taken 
from [l^]. INSET: The Knight shift versus bulk susceptibility. 



This result suggests that there are local lattice distortions 
with a symmetry lower than cubic at the Ga impurities. 

Experimental 

In order to avoid radiological contamination of the 
NMR probe, we constructed an epoxy assembly to en- 
capsulate the NMR coil as well as the Pu sample. How- 
ever, since the Pu is subject to both self-heating as well 
as rf heating due to induced surface currents by the NMR 
pulses, it is crucial to establish a solid thermal contact to 
the sample. We therefore bored out the epoxy along the 
axis of the NMR coil, and secured frits (metallic screens) 
at the ends to allow the cold He gas to enter along the axis 
of the coil and establish a thermal contact with sample. 
The mesh of the frits is sufficiently fine to prevent any 
external contamination. This design allows us to cool the 
sample down to 2K with no significant heating problems, 
as measured by the nuclear polarization. 

The Hamiltonian of the Ga nuclei is given by: 

7^ = 7ft/ • Ho + ^(3/2 - /2 + ,y(/2 _ /2)) + ^j^^p (1) 

where 7 is the gyromagnetic ratio, Hq is the external 
field, V, = 3eQKc/20, r/ = (Ka - T4b)/Kc, Q is the 
quadrupolar moment of the Ga, and Vaa are the com- 
ponents of the EFG tensor. The hyperfine interaction is 
given by Hhyp — I • X^i ' ' where the sum is over 
the various spin contributions (conduction electrons, lo- 
cal spins). For a site with cubic symmetry, Vaa vanishes, 
so the effective Hamiltonian for the Ga nuclei in (5-Pu is 
simply: 

n = 7^(1 + K)i,HQ, (2) 



where 

K{T) = K^ + Y,A,x^{T), (3) 

i 

Kq is a temperature independent orbital shift contribu- 
tion, and Ai is the hyperfine coupling to the i*'' compo- 
nent of the electronic susceptibility, x{T)- 

In Fig. ^ we show the Ga spectrum at 4K. Note 
that this spectrum consists of all three nuclear transi- 
tions {h = +i ^ = +5 ^ -h^- '^-\'^ -§)' 

since the quadrupolar Hamiltonian does not lift the de- 
generacy. The spectrum is broadened, however, by a dis- 
tribution of EFG's, as we discuss below. The center of 
the resonance is somewhat lower in frequency than that of 
pure Ga metal, indicating a different Knight shift. We fit 
the spectra between 4K and lOOK to a Gaussian; K and 
the rms second moment, tr, are shown in Figs. (|2I3() . If 
we assume a single component of magnetic susceptibility, 
then following Eq. (PJ, we extract and A by plotting 
K versus x, where T is an implicit parameter (see inset 
of Fig. ©■ We find = 0.65% and A = SOkOe/^s. 

Spin Lattice Relaxation 

The spin lattice relaxation measures the time scale for 
the nuclear spin system to acquire an equilibrium Boltz- 
mann distribution among the energy levels. Transitions 
between the nuclear levels are induced by time dependent 
fiuctuating fields that are perpendicular to Ho. The spin 
lattice relaxation rate, T^^ is give by: 

^ = ^i-oE^i(q)^^^^ (4) 

q 




FIG. 3: The rms second moment of the spectral hnewidth 
versus temperature of the ®^Ga hue in 5-Pu. INSET: The 
echo integral versus r^, where r is the pulse spacing between 
the 90° and 180° pulses. 



where A± (q) is the spatial Fourier transform of 
the hyperfine coupling in the perpendicular direction, 
X_L"(q, w) is the dynamical q-dependent susceptibility, 
and the sum is over the Brillouin zone. In this case, we 
expect that A is isotropic and q- independent, thus the 
spin lattice relaxation rate is sensitive to fluctuations for 
all q, and any magnetic correlations or magnetic order 
should be reflected in T^^- 

We measured the ^^Ga T^^ by inversion recovery. In 
Fig. Q we show data at 4K, which shows a single expo- 
nential. In principle, the recovery of the / = § Ga can be 
more complex. However, since the satellites overlay the 
central transition (strictly speaking, the Hi pulses have a 
bandwidth much greater than the EFG distribution) all 
nuclear levels are inverted by the pulses, and the recovery 
function simplifies to a single exponential of the form: 



M(i) = Mo(l-/e-*/^^), 



(5) 



where Mq, /, and Ti are fit parameters. If the EFG were 
actually much greater, and the satellite transitions {Iz = 
±1 ^ ±^) were either too broad to observe, or missing 
from the spectra for some reason, then the recovery of 
the observed central transition (/^ = +^ ^ — i) would 
be of the form: 



M(t) = Mo 



1-/ 
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(6) 



Clearly, as observed in Fig. Q), the data indicate a single 
exponential, which supports the conclusion that the Ga 
is located at cubic site. In the inset, we compare the 



FIG. 4: The magnetization decay of the Ga line at 4K, and 
fits to a single exponential decay (solid line) and the decay 
function for the central transition (dotted line). INSET: The 
magnetization recovery data for both isotopes of Ga. The 
time scale for the ^^Ga has been scaled by C^'y/''^'y)'^ , as 
expected for magnetic relaxation. 



relaxation of the ^^Ga with that of the ^^Ga. The time 
scale for the ^^Ga has been scaled by C^l/^^l)'^ (^^^ 
Eq. 10}); the fact that the data for both sets of isotopes 
fall onto the same line indicates a magnetic relaxation 
mechanism, rather than a quadrupolar one (as might be 
expected for structural fluctuations, for example). 

In Fig. (O we show the temperature dependence 
of T]"^. In conductors where the dominant relaxation 
mechanism is by scattering with the conduction elec- 
trons via a single contact interaction, such as in the al- 
kali earth metals, the Knight shift and the T^^ are re- 
lated TiTK^ — K, where the Korringa constant is 
given by: n = fj.^/iTkBfi'y^. For convenience, we define 
S = TiTK'^/k, so for the simple case, S{T) is unity. 
Clearly this quantity approaches a constant value at low 
temperatures, as seen in the inset of Fig. jSjl. This result 
suggests the absence of magnetic correlations or magnetic 
order of the Pu spins, which would likely contribute a 
strong temperature dependence to TiT. 

The fact that S{T 0) ~ 0.05 suggests that hyper- 
fine interaction in S-Pu is more complex than in the alkali 
earth metals. Electron-electron correlations can give rise 
to an increase in S{T), however 5' > 1 is seen only in 
systems with a single hyperfine coupling 0, 0. The 
fact that S{T) < 1 suggests that there is more than one 
hyperfine coupling mechanism to the conduction elec- 
trons. In particular, if there were two hyperfine cou- 
plings, Ai and A2, to two different conduction bands, 
then the Knight shift would be given by the sum A1+A2, 
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FIG. 5: The spin lattice relaxation rate versus temperature for 
the Ga (•), and Al (□, from [13]) in S-Pu. INSET: TiTK^/n 
versus T, where k is given by the Korringa constant. The 
dotted line is a fit to T^TK^/k = So{l + (T/Tof), where 
5o = 0.046, and To = 25K. 



whereas the spin lattice relaxation rate would be given 



by the sum of the squares 
might expect TiTK^ - ■ ^^i±^ 



h ^2 , so that naively one 
■ 2 , .a . If Al and A2 have differ- 
ent signs, then the measured ratio can be less than unity, 
as we observe. In fact, the measured S{T ^ 0) suggests 
A2/A1 — —0.9. A second, negative hyperfine coupling via 
core polarization is common in transition metals such as 
Platinum 14]. Given the complex electronic structure of 
(5-Pu, and the result that the f electrons have itinerant 
behavior, such a result is not surprising [isllr^. 

As seen in the inset, S{T) has a quadratic temperature 
dependence; we fit the data to S(T) = Soil + (T/Tq)^). 
In transition metals, Tq is a measure of the energy scale 
on which the density of states has structure away from 
the Fermi energy. Specifically, Yafet and Jaccarino find 
nkBTo = {p/'id?p/dE'^)Ep, where p is the density of 
states . Comparison with photoemission data or band 
structure calculations is difficult, however, since kpTp = 
2.2meV is less than their typical energy resolution |j,ll5j- 



Inhomogeneous Broadening 

We now turn to the width of the NMR line. In princi- 
ple, the resonance can be broadened by inhomogeneitics 
in iJoi by magnetic interactions with the other nuclei in 
the system, by local magnetic moments on the Pu sites, 
and by quadrupolar effects. The magnet has a homogene- 
ity of lOppm, so field inhomogeneitics can be neglected. 
Let us define cr^ as the second moment of the NMR line 
spectrum. For a nucleus with J > ^ we expect quite 
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FIG. 6: The hole-burning spectrum of ''^Ga at 4K. The solid 
line is the spectrum after a narrow inversion pulse, and the 
dashed line is the spectrum with no preliminary pulse. 



generally: 



2 



(7) 



where cr^ag the contribution from magnetic moments in 
the system, is the contribution from non-zero EFG's 
at the nuclear site, XnHo is the contribution from the 
anisotropic susceptibility and demagnetization fields in 
the powder, and cr^^, is the second moment of the nu- 
clear dipole interaction. Note that rj = ij' corresponds to 
like spin coupling (for example ^^Ga-^^Ga coupling) and 
rj ^ rj' corresponds to unlike spin coupling (^^Ga-''^Ga 
and 69Ga-239Pu) 13. These quantities can be calcu- 
lated for the fee Pu lattice with randomly located Ga 
atoms, and depend on the orientation of the field with 
respect to the crystal parameters 0. For our polycrys- 
talline sample, the broadening is an average over the unit 
sphere; since a^^, varies by less than 50% over the unit 
sphere we have chosen the (100) field direction for con- 



creteness. We find cr|g gg 



and cr|9^239 



0.009G2, cr|g 71 0.004G2, 
2.5G^. As seen in Fig. |j3Jl, the spectral 
broadening is much greater than the dipolar coupling. 
In order to determine the field dependent contribution, 
we measured the linewidth as a function of applied field, 
as seen in Fig. Fitting the data to Eq. iQ we 

find a{HQ = 0) = 4.6G, a value still greater than than 
dipolar second moment, so the NMR spectrum must be 
inhomogeneously broadened by ~ 4.3G ~ 4.4kHz. 

To test for inhomogeneous broadening we performed 
a hole burning test by saturating a narrow fraction of 
the line. We applied a low power pulse to the system 
prior to a broad-band echo sequence; the spectrum is 
shown in Fig. 0. Clearly, the full spectrum consists of 
the superposition of several narrow intrinsic lines. Note 
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FIG. 7: The rms second moment of the spectral line versus 
applied field. The solid line is a fit to Eq. (|7|l, and the dotted 
line is the estimated homogeneous linewidth. 



that the rms second moment for the "hole" is 6.7G, a 
value on the order of the field independent broadening, 
but still greater than the calculated value. However, the 
excitation pulse in this case was still greater than the hole 
linewidth (but narrower than tha full spectral linewidth) , 
so it is doubtful that the measured hole linewidth reflects 
the intrinsic linewidth. 

The inhomogeneous broadening arises either from 
magnetic moments in the system, or from a distribution 
of EFGs. In the latter case, we utilize the measured hy- 
perfine coupling to estimate the size of the putative mag- 
netic moment as 5 x 10~^/iB. Such moments would not 
be detected in bulk susceptibility measurements. How- 
ever, neither a nor K show any significant temperature 
dependence, as one might expect if these local moments 
were present in the system. 

In fact, the most likely source of inhomogeneous broad- 
ening is from a distribution of EFGs which contribute 
the the quadrupolar linewidth. One way to test for this 
is to compare the linewidth of the ^^Ga to that of the 
^^Ga. However, we find the surprising result that the 
^^Ga width is 76% larger than that of the ^^Ga, whereas 
the quadrupolar moment of the ^-'^Ga is 40% smallerl 
(Note, however, that the measured Xn's (see Eq. [TJ of 
the two isotopes do follow the ratio of the gyromagnetic 
ratios, as expected.) The reason for this discrepancy may 
be related to the fact that the sample had been held at 
temperatures T < lOOK for several days prior to the 
measurements of the ''^Ga. At these temperatures, the 
damage to the lattice inflicted by the radioactive decay 
of the Pu atoms is not annealed out, and is reflected in 
resistivity measurements . One expects that the EFG 
will reflect lattice damage, thus it is reasonable that the 
inhomogeneous broadening is due to lattice distortions. 



In this case, the 4kHz broadening would correspond to a 
distribution of e^^QVcc/h - 27kHz in the local EFG at 
the Ga site. 

Any perturbation from cubic symmetry at the Ga site 
will give rise to a flnite EFG and consequently contribute 
to the line breadth. XAFS studies of the Pu-Ga distances 
in (5-Pu indicate that the lattice contracts slightly around 
the Ga impurities |^, however it is not clear how the 
EFG will be modifled as a result. Further studies are 
needed to determine what fraction of the quadrupolar 
line broadening arises from radioactive damage and what 
arises from lattice contraction around the Ga impurities. 



Decay of the Echo Envelope 

The decay of the echo envelope can provide information 
about the like-spin coupling and the intrinsic linewidth. 
The spin echo is acquired by applying the pulse sequence 
90° — T — 180°, and occurs at a time 2r after the first 
pulse. The integral of the spin echo is plotted versus 
in the inset of Fig. Q . In the limit where cr^^ is highly 
anisotropic, so that the like-spin nuclear coupling is much 
larger in a particular direction, the form of the echo decay 
can be solved exactly in terms of (t,*^,^ 2^ ■ ^^i^ 
case, the echo decay is Gaussian: 



M (t) = Mo exp 



(2r)2 



(8) 



where T^, 



' ^ ri 



— 7 u^^. However, for the dipolar couplings 
discussed here, we are not in such a limit, and one cannot 
write down an exact form for the echo decay Nev- 
ertheless, we find that the echo decay is indeed Gaussian 
with a time constant T2G ~ 480/Lts. This value corre- 
sponds to a second moment of 0.105G^, a factor 10 times 
larger than calculated. Even though we do not have an 
exact form for the echo decay for dipolar couplings, one 
would expect a priori that the measured echo decay con- 
stant would be within an a factor of 2-3 of the like-spin 
second moment. Therefore, the fact that we find such a 
disparity is surprising. If the like-spin coupling is in fact 
larger than we estimated, then this result suggests one of 
three causes: (i) there exist indirect couplings between 
the Ga (unlikely due to their large spatial distances), or 
(ii) the effective Ga-Ga distance is smaller, as might be 
expected for Ga clustering, or (iii) the unlike Pu spins are 
fluctuating quickly, so that there is a fluctuating field at 
the Ga site that contributes to the dephasing of the Ga 
spins H^. XAFS studies of the Ga distribution suggest 
that the Ga is distributed uniformly, so the most likely 
explanation for the enhanced echo decay rate is (iii) ja]- 
In fact, we do expect that the Pu nuclear spins have a 
very fast spin lattice relaxation rate, since they must have 
a large hyperfine coupling. 



Conclusions 

Ga NMR in S-Pu provides information both about the 
local structure and distribution of the Ga atoms, as well 
as the electronic spin fluctuations. We find that for 1.7 
atomic percent doping, S-Pvl shows little evidence for 
local magnetic moments at the Pu sites, but that the 
hyperfine coupling between the Ga and the conduction 
electrons probably contains a contact as well as a core 
polarization term. Furthermore, the NMR spectrum is 
inhomogeneously broadened by a distribution of EFG's 
at the cubic symmetric Ga site. More detailed studies of 
the linewidth as a function of doping should yield impor- 
tant information about the Ga distribution. It is worth 
noting that these experiments were conducted at tem- 
peratures lower than the proposed Kondo temperature 
of 200-300K [13. Therefore, further studies at higher 
temperatures may shed light on the presence of spin fluc- 
tuations. 

This work was performed at Los Alamos National Lab- 
oratory under the auspices of the US Department of En- 
ergy. We thank J. Sarrao, J. Thompson, J. Wills, M. 
Fluss, and C. P. Slichter for enlightening discussions. 
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